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Abstract 

A critical analysis of the conformal approach to the theory of 2d 
turbulence is delivered. It is shown, in particular, that conformal 
minimal models cannot give a general turbulent solution, which should 
provide for constant fluxes of all vorticity integrals of motion. 

The problem of small-scale spectrum of two-dimensional turbulence is a pe- 
culiar problem among the variety of turbulent systems. The point is that 
dimensional considerations do not give a steady spectrum that corresponds 
to the vorticity cascade. The spectrum obtained from dimensional analysis 
is E(k) oc k~ 3 [1] which yields a logarithmic infrared divergence after substi- 
tution into the equations for the correlation functions. Kraichnan's attempt 
to save the spectrum from nonlocality by introducing the slow factor In -1 / 3 k 
attains convergence only in the first order of perturbation theory [2] while 
the next orders reveal divergences with higher powers of the logarithm: In 2 
etc. The fact that the powers of the logarithm increase with the order of 
perturbation theory suggests that a substantial renormalization of the index 
occurs. No successful attempts to work out the divergences or to show that 
they are cancelled are known to us. The existence of alternative predictions 
for the steady spectrum E(k) oc k~ 4 by Saffman [3] and E(k) oc k' 11 ' 3 by 
Moffatt [4], show that this is still an open problem. 

A fairly new approach to the problem has recently been introduced by 
Polyakov [5]. He suggested to borrow a set of correlation functions from 
conformal field theory to satisfy a chain of equations following from Euler's 
equation. Conformal invariance is assumed for the set of the simultaneous 
correlation functions (not for the Euler equation itself). This is based on 
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the idea that nonlinear interaction renormalizes all the properties in such a 
way that at criticality (for an infinite Reynolds number, in our case) system 
possesses the highest possible symmetry (absent in the initial equation). It is 
convenient to start from the Navier- Stokes equation written for the vorticity 
field u;(x, t) 

to + e a pd a ipdi3^.ip = vAoj , (1) 

ip is a stream function giving the velocity field as v a = e^dpip. Here and 
below we use the notation d/dx a = d a . 

Our aim is to find a stationary set of equal time correlation functions 
J n (xi, . . . , x n ) = (a>(xi, t) . . . u(x n , £)). The brackets denote an average with 
some time independent probability distribution: 

n 

Mxi, t) . . . £j(Xp, t) ... cj(x n , t)) = . 

p=l 

Such a stationary set is expected to exist in the inertial interval of scales, i.e., 
for distances that are much less than the scale L of an external pump and 
much larger than the viscous scale a. It is possible, then, to neglect viscosity 
in (1) using instead a careful point splitting procedure [5] for the nonlinear 
term 

e af3 d a ijj(x)dpAijj(x) =lime af3 d a ip(x + a)9 /3 A^(x - a) 

a— >0 

where "lim" implies angle averaging. To calculate different-point pair corre- 
lators like ip(x + a)^>(x — a), the fusion rule of the type 

M [if] = [0] + • • • (2) 

should be used. Here we follow the notations of Ref.5 so that [-0] means the 
conformal class of ip, i.e. itself together with the operators L_ ni . . .L^ nk ip, 
L_ n being Virasoro generators [6]. Both ip and <j) are presumed to be taken 
from a set (primary fields) of some conformal field theory. 

Polyakov suggested to use the so called minimal models [6] which contain 
a finite number of primary fields. The main subject of this paper is to show 
that a minimal model cannot serve as a general turbulent solution. 

The primary field (f) provides the main contribution in the operator prod- 
uct expansion (OPE) @ in the small-scale region, i.e. it has the smallest 
conformal dimension. The important thing is that the scaling indices (dimen- 
sions) of the fields are not additive so generally A^, ^ 2A^, and a dimension 
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defect appears. The energy density in the wave number space is expressed 
via \ip\ 2 and is 

E(k) oc k^ +1 , (3) 

To choose an appropriate solution from the wealth of conformal solutions, 
one should impose some additional conditions that follow from the symme- 
tries or conservation laws specific for the problem in question. According 
to Fjortoft's theorem (see e.g. [7]), the vorticity is the relevant quantity 
in the problem of small-scale turbulence (while the energy flux determines 
large-scale turbulence). Following Kraichnan [1] who developed a simple and 
efficient (though uncontrollable, of course) closure in terms of double corre- 
lation functions, the enstrophy 



H 2 = J to 2 (x)d 2 x 



which is a motion integral of Euler's equation is usually taken into account. 
A steady turbulence spectrum in the small-scale region should provide for a 
constant enstrophy flux over the scales which yields [5] 

(Co(x + r)uj{x)) oc r° = const . (4) 

Puting to = Aip and Co oc (L_2^-i — L^L^^cj), Polyakov obtained [5] 

(A^ + 2) + (A^ + l) = . (5) 

As one can see, the enstrophy flux is expressed through the triple correlation 
function which can be expressed by the fusion rule (0) in terms of the double 
correlation function. 

Equation (|5]) can be obtained also by requiring the rate of the enstrophy 
dissipation to remain constant while the viscosity v goes to zero: 

dH 2 v fV«.,„ r a \ . 



dt 



■ / k 2 H k dk oc ua- 6 ~ iA ^ oc v A ^ A v- . (6) 



The last estimate was given by using the expression for the viscous scale v oc 
* that follows from the comparison of the nonlinear and the viscous 



2A^-2A 



terms in the Navier-Stokes equation. The natural physical assumption that 
the main dissipation stems from small scales is valid if > —3/2 or A^, > 
A^ as it can be seen from (|]||). 



4 



Note that Kraichnan's dimensional approach would correspond to addi- 
tive dimensions A^ = 2A^ giving thus A^ = — 1 and E(k) oc k~ 3 . However, 
we have arguments that suggest that this in not a conformal solution from 
the set of minimal models (see Appendix). 

To ensure that the conformal set of correlators is a steady solution, 
Polyakov imposed an extra condition requiring that rhs of (H) (which de- 
termines the time derivatives of the correlators) vanishes in the ultraviolet 
limit. Since ip(zi)ip(z2) = (zi — 22) *~ ^(#2) + • • • as z\ — > z 2 , the following 
inequality arises: 

A > 2A^ , (7) 

thus, using (|), A^ < — 1. 

Conditions ([5]) and (0) by no means determine a single solution. The min- 
imal model (2,21) presented by Polyakov [5] is nothing but the first example 
from an infinite family. The curious reader can find the first few hundred so- 
lutions in our preprint flOf . Polyakov's solution corresponds to the minimal 
number of primary fields (in this case 10). One could, in particular, find the 
minimal model (5,72) that gives A^ = A( 125 ) = —7/6 and E{k) oc k~ u / 3 as 
in Moffatt's spectrum. 

Usually, when speaking about a turbulent solution carrying a constant 
flux, one should check that two conditions are satisfied: i) The solution 
should be local in fc-space which means that distant scales do not interact 
substantially; this should be provided by the convergence of the integral de- 
termining the flux in fc-space; ii) The constant that arise in this (converging) 
integral should be nonzero and have the correct sign to satisfy the boundary 
conditions in /c-space, i.e., the pumping and damping. Any solution of (^) 
and (^) violates both of these conditions. 

Inequality (0) means that A^ < —1 i.e. the spectrum (|) is steeper than 
Kraichnan's one (which yields an infrared logarithmic divergence) so that 
a power infrared divergence arises for any solution. To show this, let us 
consider the Euler equation in the momentum representation: 

^ = J r 12 uiu*5(k + ki + k 2 ) d 2 hd 2 k 2 . 

Here the vertex Ti2 = [kik2](A;f 2 — k^ 2 ). The equation for the enstrophy 
density h 2 (k,t)5(k + k') = (oj^k 1 ) by virtue of the conservation law can be 
written as a continuity equation h 2 + div P 2 = 0. The flux is constant in a 
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steady state: 



k 

P 2 = 2tiJ kdk J r 12 (u k u 1 u 2 }6{li + k x + k 2 ) d 2 hd 2 k 2 . (8) 
o 

This is the Fourier transform of (|]) . Infrared convergence of (§) is determined 
by the asymptotics of the triple correlation function at k\ <C k ~ k 2 . Such 
an asymptotics is expressed via the second-order correlation function ||: 

(u k uiu 2 ) = Jkuh^h^ityih/fySik + ki + k 2 ) . (9) 

Here <&(k) is a power function which value is irrelevant to the problem of 
infrared convergence. Substituting (|) into (||), one gets the following integral 

kh 2 (k)dk . (10) 

If h 2 (k) oc k~ m , then fllOD converges for y 2 < 2. It means that the energy 
spectrum E(k) is local if it is less steeper than Kraichnan's k~ 3 . For the 
conformal solutions, y 2 = — 2A^ and the convergence condition gives A^ > 
— 1. Equation @ gives the scaling exponent of the triple correlation function 
y 3 — 4 which corresponds to ([|) only in the case of convergence. 

As far as the second condition (of nonzero flux) is concerned, the three- 
point function {ipipip) ~ ((fnf>) is equal to zero since primary fields with dif- 
ferent dimensions (A^, ^ A^) are orthogonal for minimal models. The cases 
where ip appears in the operator product expansion of [ip] [ip] (like those with 
A^, = A^ = —3/2 which one can find in fllPf ) break parity and therefore 
should be excluded. One can consider models with A^ = A^ = —3/2 but 
[ip] [ip] does not contain ip in the OPE, this solution corresponds to the spec- 
trum E(k) oc k~ 5 which seems to be unphysical since the dissipation integral 
(0) is not defined by the ultraviolet region but stems from the whole /c-space 
implying thus the absence of the inertial interval in this case. 

The conformal solutions in question are thus fluxless. This is quite natural 
since they are invariant under time reflection while a nonzero flux corresponds 
to an irreversible dissipation. The second difficulty (zero flux) remedies to 
some extent the first difficulty (nonlocality) since one should not require 
the convergence of the integrals that are identically zero. Physically this 
corresponds to the fact that the spectrum of a system in thermodynamic 
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equilibrium should not be local, unlike the cascade spectrum. A fluxless 
spectrum corresponds to an equilibrium case. 

But, since we are discussing a nonequilibrium situation, there still remains 
the question, how the spectrum carries nonzero flux from the pumping region 
to the viscous region of scales. Polyakov's suggestion is that small (noncon- 
formal) deviations from the power law due to infrared cut-off (pumping scale) 
could provide nonzero values for the enstrophy flux. (It is unclear if such a 
nonlocal yet cascade solution exists until it is explicitly found by solving the 
matching problem). The physical correlators are thus assumed to be close 
to the equilibrium ones in the inertial interval of scales. This is similar to 
what happens in two-dimensional optical turbulence described by the Non- 
linear Schrodinger Equation + A\I/ + T|\I/| 2 = 0. For wave turbulence, 
the small-scale spectrum is e(k) oc k a ~ m ~ d , where d is the space dimension 
and a and m are the scaling indices of the Hamiltonian coefficients (i.e. the 
frequency and the four- wave interaction coefficient respectively) ||11|| . For the 
NSE, a — d — 2 and m = so that the turbulent spectrum is e(k) = const 
which coincides with the equilibrium equipartition. This spectrum (which is 
an exact steady solution) is fluxless too. Numerical simulations of the NSE 



show the nonequilibrium spectrum to be close to e ~ const [12]. The same 
coincidence of the equilibrium and turbulent spectra takes place for common 
turbulence of Langmuir and ion sound waves in plasma, the spectra carrying 
fluxes acquire logarithmic factors in this case [TT| . 



In our opinion, unlike the above cases of wave turbulence, conformal 
approach, namely, the assumption that the turbulent spectrum of 2d hydro- 
dynamics should be close to an equilibrium one, is not based on solid ground. 
One would like to see the degeneracy that prescribes the coincidence of the 
turbulence spectrum with an equilibrium one. 

Polyakov suggested to distort the spectra by analytical (in x-space) con- 
tributions. It corresponds to introducing 5(/c)-terms in the momentum repre- 
sentation JTB|. As a result, the flux (g) is determined by 5-terms and it is thus 
automatically independent of k. To get the condition (|[), Polyakov required 
that the flux is independent of the infrared cut off which is the pump scale. 
This condition is satisfied in the (purely theoretical) case of the excitation by 
an external force independent of the velocity field. If one consider, however, 
a more physical case of a large-scale instability (dh 2 /dt + divP2 — 7(^)^2), 
then the flux depends on the value of h 2 (k ) at the pump region. If one es- 
timates this value by extending the power solution from the inertial interval 
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h 2 (k) = P a k~ s , one has 



P oc ^- s)/(1 - a) . 



What is important here is the factor 1 — s. It has clear physical meaning: if 
one considers the integral of the total enstrophy H = J h 2 (k)dk, then 1 — s 
determines whether the enstrophy-containing region is in small or large scales. 
The enstrophy flux increases if one shifts or extends the pump toward the 
enstrophy-containing scale. Conformal solutions correspond to s > 1. There 
are thus three possibilities: 1) conformal solution cannot be matched with 
an instability increment; 2) the value h(ko) is much less than one naively 
estimated by extending the conformal spectrum; 3) h(ko) is well matched 
with a conformal part but into a pump integral you should substitute a 
much smaller value. Maybe 3) means that actually only small number of 
eddies grow due to the instability since most of large eddies are quasisteady 
and do not provide an instability leading to the direct cascade, while the 
whole set of eddies constitutes the conformal solution. 

And what is more important, the above conformal approach does not 
account for the presence of an infinite set of motion integrals 



The conservation of H n follows directly from the fact that the equation 



conserves the integral / F(uj) dxdy, where F is an arbitrary function and the 
Hamiltonian 7i is an arbitrary functional of to [not only TC — j ipudxdy, 
which gives the lhs of ([!])]. As one can see, even an infinite number of motion 
integrals does not fix the system but only its class. 

Most authors feign indifference to the existence of the infinite number of 
motion integrals in 2d turbulence. Some arguments that only quadratic inte- 
grals (i.e. energy and squared vorticity) should be taken into account while 
considering thermodynamic equilibrium were given by Kraichnan [p], |H| . 
However, an arbitrary turbulent pump generally produces a nonzero input of 
all integrals H n . The theory should describe the fate of these integrals. 

One could try to presume that H 2 cascades while other integrals jump 
from the pumping to the damping and do not influence the structure of the 




(11) 



^ + e a/3 d a —dpuj = 
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steady spectrum in the inertial interval of scales. Let us show that this is 
impossible and that the conditions for interaction locality are the same for 
any n. 

The asymptotics of high-order correlation functions can be expressed via 
h 2 (k) by the same means as asymptotics @. This can be done by using 
Wyld's diagram technique which allows one to express any correlation func- 
tion via series containing the double correlation function and Green's func- 
tion. For k — > 0, the double correlation function increases as k~ x ~ y2 while 
Green's function behaves as k~ x . Here x is a dynamic exponent, its value is 
irrelevant for our consideration. Since we are dealing with the case y 2 > the 
relevant terms in the series are those which contain the double correlation 
function of small wavenumbers 0. Such subsequences could be summed up 
to get exact expression for the asymptotics of any correlation function. Pro- 
ceeding in this way we get the same convergence condition y 2 > 2 for any n. 
In particular, Kraichnan's spectrum gives the expression for P n containing 
(In /cL) n_1 . We should thus conclude that different vorticity integrals behave 
in the same way in fc-space. 

The conservation of H n means that the Euler equation can be also repre- 
sented in the form dh n (k) / dt + div P n = so that the steady-state condition 
takes the form P n = const. Here the spectral density h n is expressed via the 
n-th correlator while the flux P n via the n + 1-th correlator. Condition of 
the flux constancy gives the exponent of the n-th correlator y n = 2n — 2 un- 
der the same condition of interaction locality as Note that Kraichnan's 
spectrum formally satisfies these conditions (formally but not really because 
of the divergences). Indeed, the dimensions are additive in this case so that 
estimating the rate of the viscous dissipation of H n (n > 2) similarly to (|6|), 
one gets 

1 L Cl (n-l)(A v , + l)+A + 2 

H n = v k 2 h n (k) dkocv A ^" A * . (12) 

Kraichnan's spectrum gives 2A^, = A^ = —2 so that fll2|) is satisfied for 
any n. Indeed, this corresponds to ip(r) oc |r| 2 so that the vorticity uo has 
zero scaling dimension (maybe logarithmic). All powers of the vorticity can 
thus have constant fluxes in /c-space simultaneously. Actually, Kraichnan's 
spectrum would equally respect all conservation laws if it was local. 

Note that we estimated the n-th correlation function in (O) in the sim- 
plest way just to show that A^, = —1 formally satisfies all conservation laws. 
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As it is shown in the Appendix, a minimal model cannot have = 2A^ = 
—2. Any conformal model that satisfies Polyakov's conditions (|^0) gives 
some negative defect of dimensions while calculating correlation functions. 
It means that the exponent in (|T2| ) is an upper bound so that by virtue of 
(H) one can get the inequality 

with some positive constant C. Therefore, for Polyakov's solution with A^ < 
— 1, the integrals H n with n > 2 are not dissipated in the inviscid limit. 
Consider, for example, the flux of H 3 : {ujuj 2 ). It is expressed via the correlator 
((fnpip). If we fuse and ip and take the most relevant operator, say, x so 
that the condition of the flux constancy is A x + 3 + + 1 = 0. For example, 
for Polyakov's solution (2,21), A x = —15/7 so that the condition of the flux 
constancy is not satisfied and dH 3 /dt — > as v — > 0. 

For the general case we need to calculate (ujuj 71 ' 1 ) which is given by 
((pip n ~ l ) . If we take ip n ~ l to be the most relevant operator, say Xn-i, m 
the OPE of [if)] . . . [tfj] then we get the condition for constant flux of H n : 

n-l 

A^ + 2 + (n — 1) + A n _i = 0, where A n is the conformal dimension of \ n . 
One can see that to satisfy these conditions for all n one needs infinite number 
of primary fields with negative dimension. 



It was first pointed out in our preprint |L(J that a minimal model cannot 



describe the spectrum with all vorticity fluxes. This was then recognized by 



Polyakov [K| who suggested that H n can flow upscales for n > 2 and down- 
scales for n — 2. This seems unlikely. A general matching problem actually 
includes two integrals of motion: energy J v 2 dr and arbitrary functional of 
vorticity / F(u)dr. The conservation of the latter follows from Kelvin's the- 
orem on conservation of circulation. The integrals H n are the Taylor coef- 
ficients of F(u). It is natural to think that energy conservation determines 
large-scale while Kelvin's theorem determines small-scale turbulence. How 
one could divide Kelvin's theorem into two parts? Higher vorticity integrals 
cannot cascade upscales since their fluxes are incompatible with the inverse 
energy flux according to Fjortoft's theorem. Indeed, to absorb a small value 
of a large-scale vorticity one should simultaneously absorb a large value of 
the velocity (that is of the energy). Finite energy flux into the small scales 
means zero vorticity flux there. 
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Another possibility suggested in [13[] is to break the symmetry assuming 
nonzero mean values of some fields. For fields with negative dimensions this 
strongly modifies the correlation functions flIE| . In particular, the condition 
P n =const acquires the form A^+2+(n—l+A n _ 1 )—A ri+1 = 0. This condition 
gives A n « n 2 /4 at n ^> 1 which is incompatible with the assumption on 
negative dimensions. 

Let us summarize. We have the set of conformal models which are at best 
some particular solutions under an exotic condition of excitation. Kraich- 
nan's spectrum formally respects all conservation laws but it is not a solution. 
How the true solution should look like? There are two important points: the 
infrared divergences and the presence of an infinite number of integrals of 
motion. Unlike naive expectation, a true spectrum should drop with k faster 
than Kraichnan's one to compensate a nonlocal gain of the integrals by the 
modes in the inertial interval directly from the pump (cf. with [[H]]). Power 
divergence that formally arises can be considered in a spirit of a field theory 
approach |18|] and be included in the flux value as well as ultraviolet diver- 
gences in the theory of phase transitions are included in the true value of 
the transition temperature. It means that we should find a small-scale tur- 
bulent spectrum under the condition of the fixed value at pump scale. This 
value just determines the flux. Nonlocality is thus inevitable in any consis- 
tent approach to small-scale 2d turbulence. Indeed, it is easy to show that 
local vorticity cascade is impossible: if cascade was local, one can show that 
the exponent of n-th correlator would be y n = 2n — 2 which immediately 
gives a logarithmic divergence. The positiveness of the anomalous dimen- 
sions (which are the differences between true y n and 2n-2) follows also from 
the inequality (see, e.g. [|19|]) | v (/) — v(0) \< cllnl providing the absence of 



singularities in 2d velocity field. Large vortices play an important role in the 
direct vorticity cascade (see |7[] and references therein). 

Kraichnan's spectrum together with higher correlation functions having 
the exponents y n = 2n — 2 seems to be a good bare solution which should 
be subjected to the procedure of an ultraviolet renormalization. Renormal- 
ization to be done should give the anomalous dimensions linearly depending 
on the order of the correlation function to provide for constant fluxes of all 
vorticity integrals [|T^] . 
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A 

Here we summarize all the conditions set by Polyakov on the solutions of a 
turbulence problem. We describe the algorithm we used in [[UJ to compute 
these solutions and give arguments why we think = — 1 is not a minimal 
model solution. 

The conditions are: 

1. A^ + A (f) = -3 

2. < -1 

3. must be the operator with the smallest dimension in the OPE of [ip] 
and [ip]. 

The minimal models are characterized by two positive co-prime integers 
(j>, q) . For each minimal model (p, q) there is a set of ^ p ~ 1 ^ g ~ 1 ^ primary fields 
parameterized by two integers (n, m), 1 < n < p — 1, 1 < m < q — 1. The 
spectrum of conformal dimensions is given by 

(nq - mpf - (q - pf 

M = 4^ ' ( } 

From this formula we see that A( njOT ) = A( p _ nj(? _ m ) and they both corre- 
spond to the same primary field. The scaling dimension of the primary field 
(n, m) is 2A( n m ). From the first condition, using (|13"D, we get - in terms of 

p _ n^m^ + n^m^ - 8 + k 
q m\ + m\ - 2 

where k is an integer defined by 

k 2 = {n^m^ + n^m^ - 8) 2 - (n^ + - 2){m 2 b + m\ - 2). (14) 

This defines both p and q since they are co-prime. From the selection rules 
for the OPE we have n^m^ odd numbers and satisfy 1 < < 2n^ — 1, 
1 < ™><f> < 2m^ — 1. Requiring A^ be minimal we get {n^q — m^p) « 0, this 
sets m^ to be the nearest odd to 



n^m^ — 8 ± J ~ ^ n 4> m i(> + + 60 

* - 2 . (15) 
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An algorithm to calculate these turbulent solutions can be as follows. Given 
n^,rrif, one can calculate from flT5|) and if k from (14) is an integer one 
can also find p and q. Using this algorithm we got a set of solutions up to 
= 10 0. 



We also want to check if there exists a minimal model for which = — 1 
or A^ = —2. We set [n^q — m^p] to its minimal value and since p, q are co- 
prime there exist n^m^ such that n^q — m^p = 1 and we get, using ([13|), 
—2 = — , H ' or 

Apq 

p 2 - lOpq + q 2 - 1 = 0. (16) 
for large p, q one can neglect 1 and solve the quadratic equation to get 

g) ± = 5±V24, (17) 

the two solutions correspond to the symmetry between p and q. (This solu- 
tion is irrational so in practice one can take p, q co-prime such that - is close 
to this value, one can get any accuracy by taking large p, q). From A^ = — 1 
and (0) we have in terms of as the nearest integer to 



= n^- + 2J-. (18) 




The sequence defined by 

a k+ i = I0a k - a k -i a = 0, a x = 1 

serves as a solution to ( |i~6l) taking p = a k , q = a k+ i for any k — 2, 3, the 
k-th term is given by 

(5 + v / 24) fe - (5 - y/2A) k 

ak ~ 27H 

Demanding that this (p, q) model has A^ = — 1 leads to the condition that 
y/4pq + 1 is an integer or, in terms of this sequence , that 

(5 + y^i) 2 ^ 1 + (5 - V2A) 2k ~ l + 14 
24 

is an integer squared, which is unlikely. 

In our preprint |TD[ , one can find the models with (p, q) close to ( |TTD and 



n^m^) satisfying ([18]): (13,129), (39,389), (109,1082), (232,2295) and 



(69, 686) with -A, = {§, 1, Hi, gfi, H§ respectively. 
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